Boundary Layer Analysis of the Ridge Singularity in a Thin Plate. 
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Large deformations of thin elastic plates and shells present a formidable problem in continuum 
mechanics which is generally intractable except by numerical methods. Conventional approaches 
break down in the limit of small plate thickness due to appearance of discontinuities in the solution 
which require boundary layer treatment. We examine a simple case of a plate bent by forces exerted 
along its boundary so as to create a sharp crease in the limit of infinitely small thickness. We find a 
separable boundary layer solution of the von Karman plate equations which is valid along the ridge 
line. We confirm a scaling argument [Q] that the ridge possesses a characteristic radius of curvature 
R given by the thickness of the sheet h and the length of the ridge X viz. R ~ /i^/^X^/^ The 
elastic energy of the ridge scales as i? ~ K{X/hY^^ where k is the bending modulus of the sheet. We 
determine the dependence of these quantities on the dihedral angle of the ridge tt — 2q. For all angles 
R ~ q"*^^ and E ~ a'^^ . The framework developed in this paper is suitable for determination of 
other properties of ridges such as their interaction or behavior under various types of loading. We 
expect these results to have broad importance in describing forced crumpling of thin sheets. 



I. INTRODUCTION 



Mechanical properties of a thin elastic sheet undergo- 
ing large distortions is a rich problem which is still largely 
unresolved despite great effort The postbuckling be- 
havior of thin- walled structures is of importance in safety 
design and in development of energy absorbing structures 
[^,|| . Some microscopic solid-like membranes can also be 
found in a crumpled state. Examples of such membranes 
include phospholipid bilayers below the 2D freezing point 

, networks of tropomyosinj ^ , and monomolecular lay- 
ers such as graphite oxide |7| or molybdenum disulfide 
[Q. In an entirely different context, formation of moun- 
tain ranges is a result of strong deformation of earth's 
tectonic plates which to an extent behave as elastic sheets 

Several types of distortion in elastic membranes have 
been recently analyzed. Thermal fluctuations of a flat 
membrane roughen the surface, but preserve its overall 
flat shape The surface is thought to lose its flat- 

ness for sufficiently high temperature, at the "crumpling 
transition." By contrast, certain defects deliberately in- 
troduced into fiat membrane destroy its flatness. The en- 
ergy associated with these distortions is well understood 
[0|. Different distortions occur when a flat membrane 
is collapsed by an external force. The incipient deforma- 
tion in response to a load and the associated buckling 
instabilities are classic subjects of continuum mechanics 

Large deformations of thin plates and shells can be 
described by two coupled quadratic fourth order partial 
differential equations whose explicit solution exists only 



in a handful of cases [ 
however, that the smal 



The main complication is, 
parameter related to the thick- 
ness of the shell usually multiplies the highest deriva- 
tive or the nonlinear term in the equations This 



leads to formation of boundary layers in the asymptotic 
limit of small shell thickness. An example in which such 
boundary layer appears is a thin plate bent by torques 
applied at the boundary flrst analyzed by Kelvin, Tait 
and Friedrichs [ pO| . Other well-known boundary layer 
phenomena which include the Prandtl's boundary layer 
in a flow at large Reynolds numbers arise under sim- 
ilar conditions. 

The novel character of the boundary layer problem 
considered in this article is in its potential relevance to 
a broad range of crumpling phenomena. A number of 
empirical studies have investigated the statistical scaling 
properties of the crumpled state |17-19|. Current the- 
oretical understanding of this forced crumpling is in its 
initial stages, however. It was argued in Rcf. [2^ that the 
structure of a thin crumpled sheet can be thought of as 
a collection of flat facets bound by a network of straight 
ridges which terminate at sharp vertices. These ridges 
have the same nature as the boundary layer arising in 
a simple geometry considered in this article. The scal- 
ing properties of the boundary layer solution allow one 
to make far-reaching conclusions about the morphology 
and the elastic energy of a crumpled thin sheet. 

In the next section we derive the equations which de- 
scribe large deflections of thin elastic plates and introduce 
the boundary value problem which exhibits the novel 
boundary layer. Scaling properties of this boundary layer 
are examined in Section III. A separation of variables 
Ansatz is introduced in Section IV. The Ansatz implies 
that the local properties of the ridge solution scale with 
the distance from the boundary in a simple way. In Sec- 
tion V we derive the scaling of the separable solution 
with the dihedral angle of the ridge a. From here on we 
will refer to a as the "dihedral angle" even though it is a 
half of the difference of the dihedral angle from tt. Con- 
fusion is unlikely to result since this is the only way a 
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is used in the paper. Corroborating numerical evidence 
is presented in Section VI. Finally, the implications of 
the thickness scaling for the crumpling problem and the 
validity of the separable solution as well as future work 
are discussed in Section VII. 



II. VON KARMAN PLATE EQUATIONS 

We start with a brief review of the theory of large 
deflections of thin plates originally due to von Karman 
p2| . The equations are usually obtained within an ex- 
pansion scheme of the full three dimensional elasticity 
equations with the thickness of the plate as a small pa- 
rameter. Small strains and a linear stress-strain relation 
are usually assumed [p3| . We however prefer to derive the 
equations in a somewhat different way which is more in- 
structive while not as rigorous as the thickness expansion 
method. 

When the thickness of the plate is small the in-plane 
stresses are much greater than the normal stresses, so 
that the dependence of the in-plane stresses on the nor- 
mal coordinate is simple |l2| . Then, all variables can be 
integrated over the thickness of the plate and it can be 
treated as a two-dimensional surface. The plate has a 
preferred material coordinate system x g D some open 
simply connected domain in R^. According to the Fun- 
damental Theorem of Surfaces [Q, to define a surface 
uniquely up to an overall translation and rotation in R"^ 
it is sufficient to specify two symmetric tensors in _D: a 
metric tensor gap and a curvature tensor Cap satisfying 
a set of relations which we will write down shortly. To 
clarify the meaning of these tensors we first note that the 
strain tensor 70,^3 is defined as the deviation of the metric 
tensor from the identity 



5q/3 — ^aP + 2"fap- 



(1) 



The sum of the eigenvalues of the strain tensor 71 + 72 
is the (2D) expansion (or compression) factor and their 
difference is the shear angle . The eigenvalues of the 
curvature tensor, on the other hand, are the inverses of 
the two principal radii of curvature of the surface. Given 
the three-dimensional position r(a:i,a;2) of the material 
point {xi,X2) the metric tensor and the curvature tensor 
are given by |2J] 



gap = (9a f) • (dpr) 
Cap = n • (dadpr), 



(2) 
(3) 



where da denotes differentiation with respect to the ma- 
terial coordinate Xa, and n is the unit normal to the 
surface. To be absolutely rigorous one must distinguish 
upper and lower indices but since we define all quanti- 
ties to first nontrivial order in the strain jap and raising 
and lowering of the indices is accomplished by applying 
the metric tensor gap, raising and lowering indices only 
affects higher order terms in jap- 



In order for the tensors gap and Cap to define a surface 
they must satisfy two relations involving the Christoffel 
symbols Tapfi which are defined in terms of the met- 
ric tensor. When the strains are small the expression is 
particularly simple: 



(4) 



The first relation which gap and Cap must satisfy in or- 
der to define a surface is the Gauss's Theorema Egregium 
p^ . It expresses the Gaussian curvature K = det Cap 
which is the determinant of the curvature tensor in terms 
of the Christoffel symbols. In the case of small and slowly 
varying strain {i.e. 9^7 ^ [d^Y) the Gauss's theorem 
reads 

K = 92r2ii — 92r2i2 + riiprp22 — ri2prp2i 

= dadpjap - V^tr 7a/3, (5) 

where = d^d^. Summation over repeated indices is 
implied. Geometrically, the Gauss's equation captures 
the intuitive notion that non-zero Gaussian curvature 
(the sheet curves in both directions) must cause the sheet 
to strain. Due to historical reasons Eq. (|^) in this context 
is usually referred to as the first von Karman equation. 

The other set of relations usually termed as Codazzi- 
Mainardi equations describes how the curvature tensor 
behaves when transported around a closed curve. Again, 
when the strains and consequently the Christoffel sym- 
bols are small these relations are simple. They say that 
if the curvature tensor is parallel transported around a 
closed curve, the change is of higher order in the strain 



d^Cap = dpC 



'pLyaf 



(6) 



These equations are a tensor analogy of the condition 
satisfied by an irrotational vector field. And just as an 
irrotational vector field can be written as a gradient of a 
scalar, there exists a scalar function / such that 



Cap — dadpf. 



(7) 



Equations Eq. (||) and (||) ensure that Cap and jap de- 
scribe a physical surface. We now consider further condi- 
tions which assure that each element of the surface is in 
mechanical equilibrium. The forces and the torques act- 
ing on an infinitesimal element of the surface Sx6y can be 
found from the the stress tensor aap and the torques Map 
Jl^ as shown in Fig. la and lb. The force tensors are 
related to the deformation tensors "fap and Cap by the 
constitutive relations for an elastic body. For sufficiently 
small strains. 
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(8) 



where Y is the Young's modulus and v is the Poisson 
ratio of the elastic material [O. Here tap is the two 
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dimensional antisymmetric tensor. There is a similar re- 
lation for the torques per unit length 



Mafj — K [Cafi + Veap€prCpr] 



(9) 



which comes from considering the stress distribution 
throughout the thickness of the plate h Here 
K = Yh? I (12(1 — v^y) is the bending rigidity of the plate 
| p2| . The in-plane force equilibrium can be written down 
by inspection (see Fig. la). It reads 



(10) 



This condition allows one to write the stresses in terms 
of a scalar potential x traditionally called the force func- 
tion: 



(11) 



This relation is a tensor analog of the divergenceless vec- 
tor field expressed as a curl of a vector potential. 
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Nomal force = 2f2-= F C 8g 



C'^) (d) 
FIG. 1. (a) In-plane stresses (per unit length) acting on 
a small element of an elastic sheet, (b) Torques per unit 
length acting on the same element, (c) One-dimensional 
derivation of the normal force on a curved element bl due 
to in-plane forces F, (d) Normal shear stresses acting over 
the sides of the elementary area of the sheet. 

The resultant normal force on the element of the sur- 
face can come from three sources. First, is the external 
distributed load which we for the time being will ignore. 
Second, the normal components of the in-plane tensions 
that act on the element are non-zero due to the curva- 
ture of the surface. It is easy to become convinced (see 
Fig. Ic for a one dimensional version of the derivation) 
that that normal force per unit area due to the in-plane 
tension is haaf3Caf3 Third, changing torques in the 
sheet produce a normal force. To see this, we define the 
normal shear stresses Qa acting over the sides of our sur- 
face element as shown on Fig. Id. The normal force on 



the element due to the variation of these shear stresses 
is clearly daQaSxSy. The normal shear stresses Qa can 
be found by balancing the moments about the x- and y- 
axes of all the forces acting on the element which gives 
Qa — dpMap- Collecting the normal forces we obtain 
the second von Karman equation 



dadpMap — ha„FiC, 



(12) 



The more familiar form of the von Karman equations 
emerges when one substitutes the potentials / and x into 
Eqs. (|) and (|l|) using Eq. (|), Eq. (|), and the defini- 
tion of the potentials. They read 



I Y 

where we have defined 



(13) 



[a, 6] = eaf_i€p„{dadf3a){dfj,d^b) 

dx"^ dy^ dy"^ dx"^ dxdy dxdy 



(14) 



Notice that [/, /] is twice the Gaussian curvature K . 

In principle, the strains ^ap and the curvatures Cap 
obtained from the von Karman equations as a function 
of the material coordinates define the surface uniquely 
(up to position in R''). The problem of finding the shape 
of the surface from the strains and the curvatures is how- 
ever highly nonlinear and intractable in general. On can 
make progress in a limited class of deformations in which 
the normal to the surface does not change much. In that 
case the so called Monge coordinates are appropriate. 
The undeformed plate is located in the x — y plane so 
that upon deformation the point originally at {xi,X2,0) 
moves to (xi -\- ui,X2 + U2,'w) in the 3-dimensional em- 
bedding space where Ua and w are functions of Xa ■ If the 
derivatives of Ua and w are small everywhere, only the 
lowest non-trivial order terms in those derivatives can be 
kept in the expressions for the strains and the curvatures. 
We obtain || 



and 



Cap = dadpW. 



(15) 



(16) 



It is clear that the normal displacement w in this case is 
precisely the potential function / defined in Eq. (|^). 

To complete the description of the deformation of the 
thin elastic plate one needs to be able to calculate the 
elastic energy stored in the sheet and specify the bound- 
ary conditions on the functions / and x- We begin by 
considering the work done on the small element of the 
surface 5x5y by the surrounding parts of the plate when 
the strain in the element changes by S^ap- That work 
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is haa/sSjapSxSy Q. The stretching energy Eg in the 
plate element is found by integrating the strain from 
to its value 7^/5 while keeping in mind that the stresses 
are proportional to the strains which introduces a factor 
of 1/2. Thus the total stretching energy in the plate is 
given by 



h f 



(17) 



One can similarly show that the work done by the torques 
Mctfi in bending a surface element is ^Ma/iCa^SxSy so 
that the total bending energy Ef, in the plate is (lil 



Eb=^ J dxdy MapCaf: 



(18) 



Expressing the strains, stresses, torques and curvatures 
in terms of the potentials x ^-nd / we get after some 
algebra 

= dxdy Md^dpx)f - 2(1 + u)dei{d^dpx) (19) 

Eb ^'^ I dxdy Md^dpf)f -2{l-v)del{d^d0f). (20) 



71 - 2a 



FIG. 2. Long strip of width X bent through a dihedral 
angle a. 

We proceed to define the boundary value problem 
which exhibits the ridge singularity. We consider a strip 
(a;,?/) S {—X/2,X/2) x (—00,00) of uniform thickness 
h made of isotropic homogeneous material with Young's 
modulus Y and Poisson ratio v. Normal forces are ap- 
plied to the edge so as to bend the strip by an angle 
n — 2a (see Fig. 2). The membrane stresses CTq^ as well 
as the torques Map vanish at the boundary (except for 
the singular point y = 0). In terms of the potentials / 
and X it means that 



dadpf = dadpx = at 



±- 



X 



(21) 



One can show that the conformation r(a;i ,2:2) of the sheet 
which is a solution of the von Karman equations mini- 
mizes the elastic energy while satisfying the boundary 
conditions |23| . We note in passing that many deriva- 
tions of the von Karman equations (see e. g. Ref. |11|) 
start with writing down these energies in an add-hoc way 
and then taking a variational derivative with respect to 
the sheet shape. 

We are left to address the question of the boundary 
conditions. As pointed out in Ref. |^^, specification of 
the boundary conditions for the von Karman equations is 
tricky and is rarely done in a completely rigorous or even 
correct manner. The main problem arises when a fixed 
shape of the boundary in the embedding space is pre- 
scribed. Since the equations are in terms of the material 
coordinates, the knowledge of the shape in the embed- 
ding space requires solving the equations first. There- 
fore, a boundary condition that makes reference to the 
embedding space is intractable. Only in the case when 
the material coordinates can be replaced with the Monge 
coordinates mentioned above can one prescribe a fixed 
shape the boundary. In other cases the best one can do 
is to specify the curvature tensor and the strain tensor at 
the boundary of the material coordinates domain. One 
can imagine pulling or pushing on the edge of the plate 
which amounts to specifying the stresses at the edges. 
Torques Mafj and normal forces Qp can also be applied 
at the edge. 



The condition that the strip is bent translates into spec- 
ifying the curvature tensor at the boundary. We assume 
that the boundary is a geodesic of the surface which is 
reasonable in the small strain limit. In the direction 
along the boundary the curvature is zero except for a 
sharp peak at the origin y = Q where the strip is bent. 
The width of the bent region must be of the order of the 
thickness of the plate and the curvature is of the order 
of the inverse thickness. It is convenient then to set the 
curvature along the boundary to a (5-function since the 
all the length scales in the problem are much larger than 
h. This leads to a particularly simple condition on / at 
the boundary 



f{±X/2,y)^a\y\ 



(22) 



up to an arbitrary linear function of Xa ■ To see that the 
coefficient a is identical to the bending angle introduced 
above consider the following integral along the boundary 
(or any other geodesic which approaches x = const) 



9f, 



00) 



dy 



dy'^ 



dy 



dh 

dy 



a, 



(23) 



where we have used the definition of the curvature tensor 
Cyy — dyh-ty 'dud thc fact that the integral is taken along 
a geodesic. Here fi is the unit normal and ty is tangent 
vector to the surface in the y-direction ||2^ . The geomet- 
ric meaning of this integral is the length of the contour 
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which is inscribed by the end of the normal vector fi on 
the unit sphere as n is transported along the geodesic. 
Using this argument and a reasonable belief that the so- 
lution far away from the ridge approaches that of an un- 
perturbed flat strip on can deduce that the derivative of 
the potential / along the transverse direction y does not 
depend on x far away from the ridge. This fact will later 
allow us to gain insight into the nature of the solution. 

The physical reason one expects non-trivial singular 
behavior in the small thickness limit is the interplay be- 
tween the bending and the stretching energies. The key 
observation is that the bending modulus k vanishes faster 
with the thickness of the plate (as h^) than the stretching 
modulus Yh [|2| so that thin plates bend in preference to 
stretching. Therefore, in the limit of zero thickness, the 
bent plate develops a sharp crease at which the solution is 
irregular (the curvature in the y-direction diverges). At 
a finite thickness, as shown in ref. [|], the crease softens 
to a characteristic width governed by the competition of 
the bending and the stretching energies in the region of 
the crease. 



III. BOUNDARY LAYER SCALING 



The dominant terms in the A limit must be of the 
same order on either side of the equations. This leads 
unambiguously to /3 < which is in agreement with the 
intuitive guess that the width of the boundary layer must 
go to zero as the thickness of the sheet vanishes. Balance 
of the dominant terms gives 

and 5 = 1 (29) 

in agreement with the scaling argument of Witten and Li. 
With the knowledge of the scaling behavior of the func- 
tions / and X can obtain the small thickness behavior 
of other quantities such as, for example, the transverse 
ridge curvature Cyy — d"^ j jdy^ and the midridge longitu- 
dinal strain = (1/F)(cra;a; - vuyy) ~ {l/Y)d'^xl dy"^ 
where the expressions are evaluated at a; = and y — 
0. The leading order behavior of the transverse stress 
Uyy = d^xl (^^^ is of higher order in A and it thus can 
be ignored in the expression for ^xx- Substituting the 
rescaled quantities we find that 



We first put the von Karman equations (13) in the 
non-dimensional form by defining 

f - X y 
— . T. — — and y = — 
X ^ X 



X 



X 



f 



to obtain 



X' 



(24) 



[xJ~] 



with the small dimensionless geometry parameter 

h\ 1 
X 



A.^ 



X 



(25) 



(26) 



multiplying the highest derivative. The solution to the 
so-called reduced problem defined by A = can be read- 
ily obtained. There is no cost for bending and thus the 
strip froms a sharp crease f{x,y) = a\y\. The properties 
of the boundary layer which corrects the discontinuity in 
the reduced solution be found by the rescaling 

A''/, x = X'x, i = A"x, y = XPy. 



f 



(27) 



Note that x remains unchanged by the rescaling trans- 
formation and / and y are rescaled by the same factor to 
satisfy the boundary condition (|2^. The rescaled equa- 
tions read: 



X2~5 



d'^xd^y 

d'x 
d'^xd^y 



d^y 
.Aisd^X 



= ^^-'[xJ] 

(28) 



Cy 



X 



dy 



r,2 



and 



= A^/3 



dy 



T l2 



(30) 



These expressions give the asymptotic behavior of Cyy 
and ^xx since the rescaled quantities do not depend on A 
in the small A limit. 

Note that the width of the boundary layer has the same 
leading order behavior as the radius of curvature of the 
plate at the center of the ridge and as the "sag" which 
is the vertical deflection of the ridge shape from that of 
perfectly sharp crease. This deflection is given by /(O, 0) 
in the small dihedral angle limit a <C 1. The sag can be 
found for a general dihedral angle a from the following 
argument. Due to the a; — > — a; symmetry of the solution 
the ridge line y = is a geodesic of the surface. There- 
fore we can relate the vertical deflection of the sheet w{x) 
along the ridge to the curvature Gxx = d^f/dx^. To low- 
est order in the strain we obtain 



dx"^ ' 



(31) 



Since Cxx vanishes in the limit of small thickness, the 
leading term in the A^^^ expansion of w"{x) scales with 
the same power of A as d^f/dx^. Thus the ridge sag 
^(O) is of the same order as the ridge width and as the 
midridge transverse radius of curvature. 

We can also find asymptotic behavior of the bending 
and stretching energies of the sheet. We first note that 
in the expressions for the energies Eqs. (foh and (120) the 
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terms involving det(9c(9^/) and det{dadgx) can be ex- 
pressed as an integral over the boundary [|l2| . These inte- 
grals vanish identically for the boundary conditions con- 
sidered here. Substituting the rescaled quantities which 
remain finite as A ^ into Eqs. (ttOh and (EQ) we obtain 



(32) 
(33) 



in agreement with the energy scaling argument of Wit- 
ten and Li Q]. For a fixed thickness, these energies grow 
qualitatively slower (as X^^^) than the energy of a sharp 
crease of size X which grows linearly with X . 



IV. A SOLUTION TO THE REDUCED 
EQUATIONS 



The solution to the rescaled equations (19) can be 
sought as an expansion in powers of A^^'^, 



/ = /o + A2/-Vi + A4/^/2 + 
X = Xo + A^/^xi + A4/3x2 + 



(34) 



Plugging the series into Eq. (g8|) and matching the coef- 
ficients of powers of A^^'^ we can obtain equations for all 
orders in the expansion Eq. (p4). We get at zeroth order, 



or 



= [xo,/o] 



-i^[h, /o] 



(35) 



We note in passing the equations for /i and xi- They 
read 



5Vi 
dy^ 



. 9Vo 

dx^dy"^ 

dx^dy^ 



[xo,,/i] + [xi,/o] 
-[/o,/i]. 



(36) 



At this point we must draw the readers' attention to the 
fact that while / and x satisfy the boundary conditions 
Eq. ( pH ) or Eq. ( p2| ) there is no reason to expect that 
/o and xq do. In other words, it is likely that the ex- 
pansion in powers of A^/"^ does not converge uniformly, 
so that for a fixed A higher orders in the expansion Eq. 
( ^ ) become increasingly important as the boundary is 
approached. A solution of the zeroth order equations 
( ^5| ) might therefore be a good approximation to / and 
X only in a restricted area around the ridge away from 
the boundaries. 

To help solve equations (^5|) we observe that the 
boundary conditions do not introduce another length 
scale to the problem. Therefore, the transverse profile 



of the ridge ought to scale with the distance from its 
midpoint i = 0. We must find a functional form of /o 
and xo such that the different x — const slices are related 
to each other by a scale factor q{x) depending on x only. 
By assuming scaling in x we hope to decouple it from the 
new transverse variable ^ ~ y/q{x). We assume that 



/o(i,0-9XO, Xo(i,0 = '7>2(0 



(37) 



As we show in the appendix, variables separate for only 
one choice of 77 = /i = 1 which means that the width of 
the ridge solution scales with the distance from the vertex 
in the same way as the ridge sag. As mentioned above, 
the boundary conditions on / require that its transverse 
derivative df /dy approach a constant independent of x. 
If a similar condition is to hold for /o then inescapably 
fi — 1. The validity of this scaling hypothesis needs to 
be corroborated by some other means. In Section VI 
we present numerical evidence supporting this scaling 
Ansatz, here we only note that this evidence is convinc- 
ing. 

Substitution of the Ansatz Eq. (|3^) into the zeroth or- 
der equations (^ gives 



pT = 
vT = 



q"q^ [PiP2 



' P2P1 



mp'2+P2P'i)] 



(38) 



Separation of variables provides the equation for the scale 
factor q"q^ = A, where A is some constant. This equa- 
tion, together with the condition that q{x) be even can 
be solved by a substitution g(q) — q'ix) to yield. 




(39) 



where the separation of variable constant is related to 
q{0) via A — — 7r^(7^(0)/2. The scale factor q{x) has a 
singular derivative near the vertices, since for x = 1/2 — e 



'Z(^) 



2/3 



2/3 



(40) 



which reinforces the suspicion that the series expansion 
( p^ ) does not converge uniformly. 

Therefore to test the predicted functional form of the 
scale factor one ought to look a local property of the 
numerical solution away from the boundary and find its 
limit as A ^ 0. We discuss such a solution in Section 
VI. Numerically one has access to only a limited range 
of thickness to size ratios. (A cannot be too small other- 
wise lattice effects become significant.) Nevertheless, the 
Taylor expansion coefficients of the scale function q{x) 
around a; = can be successfully extracted from the nu- 
merics. Since there is one adjustable parameter q{0) we 
must look at the first two coefficients in the expansion 
q{x) = g(0)(l + b2X^ + h^i^ + ...). Expanding Eq. (|3|) 
we find 62 = 7r2/4 ~ 2.47 and 64 = 7r'*/48 ~ 2.03, whereas 
the numerics gives 62 = 2.53 ± 0.04 and 64 = 2.09 ± 0.07 
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where the errors are obtained from the scatter of the 
data. 



The transverse equations (38) can be integrated once 
by defining 



MO^PpiO-^p'piO for /3=1,2 



We obtain 



1 



A(t){ + D2 



(41) 



(42) 



Where Dp are the integration constants aUowing for a 
non zero asymptote of the functions (pp. The finite hmit 
of these functions as ^ — > 00 would imply a finite stress 
and curvature far away from the ridge. However, since 
we expect the series Eq. (|^ to converge non-uniformly, 
the large ^ behavior of (pp is likely to be unphysical. The 
functions pa can then be readily found from 



pm = Gpi+^ 



MO 

e ' 



(43) 



where Gp are integration constants. A further analysis of 
the transverse equations will be done elsewhere. Here we 
only remark that the simplicity and the apparent sym- 
metry of the Eqs. (42) must point to some physical sym- 
metry in the problem which forces this behavior to occur. 
Perhaps further light can be shed onto the physical mean- 
ing of the functions 4'p by writing down the expression 
for the elastic energy of the sheet in terms of the separa- 
ble solution. The leading term in the A'^/'^-expansion of 
the total energy as in Eq. ( ^3|) reads 



dx 
q{x) 



d4>i 
di 



(44) 



We can also make contact with the boundary condi- 
tion f ^ a\y\ as y ^ 00. If the integral in Eq. (^ ) 
converges, this boundary condition gives Gi = a. Now 
one can use the y —y symmetry of the problem and 
sufficient smoothness of functions involved to conclude 
that p'l (0) = which leads to 



dC dMC) 

C dC 



(45) 



We begin by noting that there exists a two-parameter 
family of transformations which produce new solutions. 
For example if 0(^) is a two-component solution, then 
tM£) = Si(j){S20 is also a solution of the equations 
( p2| ) but with a different separation of variable constant 
A' = S2/S1A. We must allow for variation of A since 
it will undoubtedly depend on the dihedral angle. Using 
the boundary condition on 0i Eq. ( ^5|) we can find the 
corresponding dihedral angle a' = SiS2a. With this con- 
dition, there is a one-parameter family of scale transfor- 
mations which produce a solution for the dihedral angle 
a' given that for the dihedral angle a. Let us fix a ref- 
erence solution (j)Q with a = 1 so that all other solutions 
are labeled by the scale factors Si and 6*2. We can now 
find how various quantities of interest are affected by the 
scale transformation. We start with the simplest cases 



a 

Pi ^P2 



S'iiS'2 
^1 

0-1/3 c-2/3 



(46) 



The relevant quantities such as the ridge sag qpi, the 
transverse curvature in the middle of the ridge {l/q)pi 
both evaluated at i = ^ = and the elastic energy 
Eq. ( ^4[ ) all turn out to depend only on the product of 
the scaling factors 5*1 6*2 which determines unambiguously 
their scaling with the dihedral angle a. Using Eqs. ( |45| ) 
we get 



qpi 

(1/9K 

dx f d^ f d(j) 

e 



q{x) 



di 



:l/3^-2/3, 
c.1/3 c-2/3 c2 03 



a 



,2/3 

4/3 



v7/3 



(47) 



One can also predict how the width of the ridge which 
scales with S'^^ depends on the dihedral angle a using 
one extra assumption. If the transverse curvature decays 
exponentially away from the ridge then the width of the 
ridge scales as the radius of curvature times the bend 
angle a. Therefore 



ridge width ^ S2 



.-1/3 



(48) 



This equation fixes the dependence of the scale factors 
Si and S2 on a so that solutions for all dihedral angles 
can be generated from a single solution for some a. All of 
the a-scaling predictions have been convincingly verified 
by the numerics. 



V. DIHEDRAL ANGLE SCALING OF THE 
SEPARABLE SOLUTION 



VI. NUMERICAL SIMULATION OF THE 
STRETCHING RIDGE 



Even without having solved the transverse equations 
one can deduce the behavior of the solution for different 
dihedral angles a just from the form of the Eqs. (42). 



We modeled an elastic sheet as a triangular lattice of 
springs of unstretched length a and spring constant K af- 
ter Seung and Nelson . Bending rigidity is introduced 
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by assigning an energy J(l — rii -112) to every pair of adja- 
cent lattice triangles with normals fii and n2. Seung and 
Nelson showed that when strains are small and radii of 
curvature are large compared to the lattice constant, this 
model membrane bends and stretches as a sheet of thick- 
ness h — ay^8J/K made of isotropic homogeneous mate- 
rial with Young's modulus Y = 2Ka/{h\/i) and Poisson 
ratio v = 1/3. The bending modulus is k = J\/3/2. We 
use a conjugate gradient energy minimization routine to 
obtain a sequence of relaxed shapes with the varying geo- 
metric parameter A. This lattice model has been used by 
the author and others to verify the A-scaling Here 
we explore the other two types of scaling which we are 
convinced govern the ridge properties. 



1 








Cyy(y) 


\ 
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\ 
















FIG. 3. Transverse curvatures Cyy{x,y) for a; = (cir- 
cles), X = b (squares), a: = 10 (triangles) and a; = 15 
(diamonds) each scaled by Cyy{x,Gi) vs. the transverse 
coordinate y scaled by Cyy{x,Q). The curvatures are 
found numerically from a 50a by 500a strip bent by a 90 
degree angle. 

A rather suggestive demonstration of the q(a;)-scaling 
is presented on Fig. 3. We plot the transverse curva- 
tures for different slices x = const obtained from a sim- 
ulated strip of dimensions 50a by 500a bent by normal 
forces applied to the x = ±X/2 parts of the boundary to 
make a 90 degree angle. The dimensionless thickness is 
A ~ 10"'^. The curvature for each slice is scaled by its 
value at the origin y = 0. The transverse coordinate is 
scaled by the inverse of the curvature at y = 0. These 
profiles collapse onto a single scaling curve. A more 
quantitative test is based on the prediction for the varia- 
tion of the transverse radius of curvature along the ridge 
{i. e. at ^ = 0) Ryy — 1/Cyy ^ q{x) . The coefficients 62 
and &4 in the expansion Ryy{x) ~ Ryy{0){l + b2X^ + b4x'^) 
can be extracted from the numerics and compared to 
those obtained from Eq. (^. 

Here we can also test the dependence of the scaling on 
the boundary conditions which cause the ridge singular- 
ity. Instead of using a long strip bent by normal forces ap- 
plied to the boundaries we can consider any shape which 
in the limit of zero thickness exhibits a sharp crease. One 



such shape chosen due to ease of implementation is a reg- 
ular tetrahedron which was previously used to verify the 
A-scaling plf . A picture of one such minimum-energy 
tctrahcdral surface is presented in Fig. 4. Shading is 
proportional to the local stretching energy. 




FIG. 4. A tetrahedron of side X = 100a and thick- 
ness h = 0.063a. Shading is proportional to the local 
stretching energy. Note the "sagging" of the ridge. 



T 1 r 




0.0003 



X 

FIG. 5. Scale factor expansion coefficients 62 (circles) 
and bi (squares) as extracted from a least squares fit to 
the functional form of the transverse radius of curvature 
along the ridge for tetrahedra oi X = 100a. 

Figure 5 displays a plot of the coefficients 62 (circles) 
and 64 (squares) versus the dimensionless thickness A ob- 
tained by a least squares fit to the functional forms of 
Ryy in the range of \x\ < 0.2 for tetrahedra of varying 
thicknesses but fixed edge length X = 100a. We see that 
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these coefBcients have a weU defined limit as A ^ which 
is approached algebraically. The extrapolations to A = 
give &2 = 2.53 ± 0.04 and 64 = 2.09 ± 0.07 as compared 
to the prediced values of 62 ~ 2.47 and 64 ~ 2.03. 

We tested the dihedral angle scaling by making a long 
rectangular strip of dimensions 50a by 200a and applying 
normal forces to its long boundaries so as to bend it by an 
angle a. For all angles we fixed A~5xl0^*. The results 
are displayed in Figures 6 and 7. Fig. 6 is a plot of the 
total elastic energy in units of the bending modulus k vs. 
the anticipated scaling variable a'/^. To a good preci- 
sion the energy does indeed exhibit the predicted scaling 
behavior. The deviation at the small bending angles is 
due to the finite size effects. Since the width of the ridge 
diverges as a ^ larger sheets are needed for smaller 
angles to avoid finite size effects. In Fig. 7 the midridge 
curvature Cyy{Q) in units of is plotted against a^/^. 
The data agree well with the scaling prediction. 




FIG. 6. Total elastic energy in units of the bending 
modulus K of 50a by 500a strips of thickness h — 0.063a 
bent by an angle a. 

Another test relevant to the question of the bound- 
ary condition dependence of the scaling predictions is 
the comparison of the coefficients R for different bound- 
ary conditions in the asymptotic form of the total elastic 
energy of the ridge £'/«; ~ i? A^^/'^a^^^. In a preceding 
work pl| ] we found this coefficient for the ridge appearing 
in the tetrahedral shape described above. Its value, which 
was found by examining the dependence of the energy on 
A for a fixed dihedral angle, is Rtet = 1.161 ± 0.003. For 
the long strip we have found Rstrip — 1.24±0.01 by fixing 
A and varying a. Here the error range reflects only the 
uncertainty arising from scatter in the numerical data. 
However, there are additional errors resulting from cor- 
rections to the asymptotic scaling not properly accounted 
for in our crude fitting procedure. We therefore suspect 
that the scaling coefficient R depends insensitively on the 
boundary conditions and varies only by a fraction of its 
value. A stronger claim of boundary condition indepen- 
dence of R would require a more convincing demonstra- 



tion. 




FIG. 7. The midridge curvature in units of for 
the 50a by 500a strips of thickness h = 0.063a vs. the 
anticipated scaling variable a^^^ . 



VII. DISCUSSION 

In this paper we described the ridge singularity in thin 
elastic plates which arises under certain types of bound- 
ary conditions. We have found three types of scaling in 
the ridge solution in the asymptotic limit of small thick- 
ness of the elastic sheet. The solution scales with the 
thickness of the sheet in a non-trivial way. In addition, 
it might scale with the distance from the center of the 
ridge and with the dihedral angle of the ridge. To as- 
certain the importance of the results we ought to discuss 
the degree to which the nature of the ridge singularity is 
independent of the boundary conditions. First we note 
that the equations used to describe the behavior of the 
elastic sheet assumed linear stress-strain relations. Since 
the strains in the ridge were found to vanish in the limit 
of small thickness, the results should be applicable to the 
real materials in which the linear stress-strain relations 
are more accurate at small strains. 

The A-scaling which was anticipated from previous 
works m2^ has been put on a more rigorous footing. 
It seems to be inevitable and depends only on the ge- 
ometry of the singularity. It is by no means a unique 
type of scaling, since Scheidl and Troger found that 
a ring ridge which appears in strong buckling of a spher- 
ical shell has a width which scales as A^/^ as opposed to 
the \^/^ scaling for the straight ridge singularity exam- 
ined in this article. The property that makes the straight 
ridge important is that its energy grows slower with size 
than that of the ring ridge which makes the ring ridge 
unstable to breakup into straight ridges when the cost of 
distorting the rest of the sphere can be overcome by the 
energy gain from the formation of the straight ridges. We 
believe therefore that straight ridges will be the dominat- 
ing type of the singularity in crumpled thin elastic sheets. 
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The morphology of a crumpled sheet can thus be repre- 
sented as a network of straight ridges. Since in the limit 
of small thickness the elastic energy is concentrated into 
a progressively smaller fraction of the crumpled sheet the 
ridges can at least in the first approximation be treated 
as independent. Therefore, once the ridge network is 
given characterization in terms of the distribution of the 
lengths and the dihedral angles of the ridges, the elas- 
tic energy of the crumpled sheet is given by the sum 

The second type of scaling found in the ridge singu- 
larity is scaling with the distance from the center x = 0. 
By assuming this scaling we are led to a separable de- 
pendence on the longitudinal and transverse position. 
Neither the scaling nor the separability, however, have 
been proven. We have only found suggestive numerical 
evidence supporting separability by examining different 
slices of the ridge as in Fig. 3. In addition, the first 
two coefficients in the Taylor expansion of the scale fac- 
tor q{x) have been found to agree with the numerics. 
A more rigorous test of separability is still needed. The 
transverse profile of the ridge is most probably dependent 
on the boundary conditions. In Ref. the author and 
others found that the degree to which the transverse cur- 
vature oscillated away from the ridge in the tetrahedron 
shape is different from that in the strip shape considered 
in this work. The simplicity of the transverse equations 
also needs to be better understood. Usually such sim- 
plicity reflects some underlying physical symmetry in the 
problem. 

The third type of scaling exhibited by the ridge so- 
lution is the dihedral angle scaling. This scaling only 
depended on the form of the transverse equations and 
not on the nature of their solution and thus is likely to 
be independent of the boundary conditions. Numerical 
evidence in Fig. 6 and Fig. 7 supports this conclusion. In 
addition, we found the same dihedral angle scaling in the 
strip shape as in the rhombus "kite" shape of Ref. . 

In future work we plan to analyze the transverse equa- 
tions and relate them to the boundary conditions in order 
to clarify whether they possess boundary condition inde- 
pendent features. A study of how non-asymptotic effects 
affect the validity of the separable solution is also in or- 
der. The framework developed in this paper is suitable 
to study of ridges under loading which is relevant to me- 
chanical properties of crumpled sheets. One could also 
introduce thermal shape fluctuations to study how they 
affect the scaling properties of the singularity, in particu- 
lar, whether the singular behavior is robust with respect 
to introduction of the fluctuations. 
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APPENDIX: 

To show that Eqs. (js^) only separate for one choice of 
the scaling exponents ji and 77 in the scaling Ansatz Eq. 
(^) we carry out the substitution for general exponents. 
We obtain 

Hl' f<f {^(V - 1KP2 + - l)P2Pi~ 

+ (^,)2^2^-, {(^- l)2(p;)2 _^(,,_ l)p';p^} . 

(Al) 

The separation-of-variables conditions on the scale func- 
tion q{x) demand that = Ai and g"g2M-r;+i ^ 
which requires that fi ^ rj. Unless the factors in curly 
braces vanish, separability also demands that {q'Y'q^ — 
A3, and (g')^^^^"'' = A^. From the Ai and A3 condi- 
tions it follows that q" q — {Ai/A3){q'y' which implies 
that q' (X q^^/^'^. Since we require that q' {0) ~ and 
(7(0) 7^ 0, one of the following two conditions must be 
true. First Ai may vanish which leads to an unphysical 
choice q" = 0. Second, the factors in curly braces may 
vanish which happens only if /i = r; = 1. We conclude 
that indeed /i = 77 = 1 as asserted in the text. 
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